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Abstract  of  Dissertation  Presented  to  the  Graduate 
Council  of  the  University  of  Florida  in  Partial 
Fulfillment  of  the  Requirements  for  the 
Degree  of  Doctor  of  Philosophy 

ATTRACTION  OF  INSECTS  TO 
ODORANT  SOURCES  IN  A  WAREHOUSE 

By 

Richard  Wendell  Mankin 

August  1979 

Chairman:   P.  S.  Callahan 

Major  Department:   Entomology  and  Nematology 

Methods  are  presented  for  determining  the  behavior  of 
an  insect  stimulated  by  an  attractant  source  in  still  air, 
turbulent  air  of  zero  average  velocity,  and  turbulent  or 
laminar  air  currents  of  constant,  nonzero  average  velocity. 
Plodia  interpunctella  (Hubner)  (IMM)  is  the  chief  test  in- 
sect for  application  of  the  methods.   To  use  the  IMM,  it 
was  first  necessary  to  measure  the  effect  of  different  con- 
centrations cf  the  sex  pheromone  (Z^E) -9 , 12-tetradecadien- 
l-ol  acetate  on  the  upwind  anemotactic  behavior  of  the  male. 
A  bioassay  was  performed  to  obtain  stimulus-response  regres- 
sion lines  at  23°C  and  34°C.   The  regression  lines  were  ana- 
lyzed by  a  new  procedure  that  accounts  for  control  respon- 
ses in  the  absence  of  pheromone  and  also  peak  responses  well 
below  100%  at  pheromone  concentrations  considerably  above 
the  lowest  detectable  levels.   On  the  basis  of  this  analvsis 


the  upwind  anemotactic  threshold  is  1.34  x  10   molecules/cm3 
at  23  C  and  1.65  x  10   molecules/cm^  at  34°c.   Departures 
from  the  2  lines  occurred  at  the  highest  pheromone  concen- 
trations  tested,  near  10  molecules/cm  .   This  suggests 
that  the  upwind  anemotactic  behavior  changes  qualitatively 
above  an  altered-behavior  threshold  about  2  orders  of  magni- 
tude above  the  upwind  anemotactic  threshold.   The  decreased 
response  at  23°C  compared  to  34°C  suggests  that  flight  in 
response  to  pheromcnai  stimulation  is  inhibited  at  cool 
temperatures. 

Calculations  of  the  above  methods  using  the  I MM  thres- 
holds and  similar  thresholds  of  other  insects  indicate  the 
following:   (1)  In  a  warehouse  a  searching  insect  is  likely 
to  be  attracted  to  a  calling  insect  if  it  comes  within  an 
attraction  space,  a  sphere  surrounding  the  calling  insect, 
ranging  from  0.4  to  2.4  m  in  radius.   (2)  The  attraction 
spaces  of  typical  sex  phercmone  traps,  emitting  pheromone 
at  rates  greater  than  0.01  ng/sec,  extend  beyond  the  boun- 
daries of  a  10  x  10  x  10-m  warehouse.   (3)  The  searching 
behavior  of  an  attracted  insect,  is  likely  to  be  altered 
from  an  extensive  to  an  intensive  pattern  if  it  comes  within 
a  calling  insect's  altered-behavior  space,  a  sphere  6-60  cm 
in  radius.   (4)  The  altered-behavior  space  of  a  trap  emit- 
ting 0.76  ng/sec  extends  beyond  the  boundaries  of  a  10  x  10  x 
10-m  warehouse.   (5)  Pheromone  does  not  fall  unless  it  is 
emitted,  along  with  a  large  amount  of  a  high-vapor  pressure 
solvent.   The  calculations  are  used  in  support  of  the  following 


(1)  The  effect  of  an  adsorptive  surface  on  the  odorant 
concentration  after  an  extended  period  of  emission  is 
negligible  except  at  positions  near  the  surface.   (2)  Traps 
with  odorant  sources  of  small  dimensions  have  greater  trap- 
ping efficiency  than  otherwise  identical  traps  with  sources 
of  large  dimensions.   (3)  The  function  of  the  altered-be- 
havior  threshold  may  be  to  increase  the  probability  of  a 
stimulated  insect  finding  a  calling  insect.   Additional 
applications  and  hypotheses  are  also  presented  for  conditions 
outisde  a  warehouse. 


VI 


INTRODUCTION 

The  probability  of  an  insect  finding  an  attractant 
source  is  determined  by  the  pattern  and  intensity  of  its 
searching  behavior,  which  are  strongly  affected  by  the 
attractant  concentration  and  the  dynamics  of  the  airflow 
(Roelofs,  1975;  Shorey  and  McKelvey,  1977) .   Several  dif- 
ferent mathematical  expressions  of  these  relationships 
have  been  presented  in  insect  attractance  models  by  Wright 
(1958),  Bossert  and  Wilson  (1963),  Bossert  (1968),  Hart- 
stack  et  al.  (1976),  Kirooka  and  Suwanai  (1976),  Aylor  (1976), 
Aylor  et  al.  (1976),  Nakamura  (1976),  Nakamura  and  Kawasaki 
(1977),  and  Roelofs  (1978).   None  of  these  models  considers 
jointly  3  problems  frequently  encountered  in  a  warehouse  en- 
vironment:  deposition  of  attractant  onto  exposed  surfaces, 
restricted  dispersal  of  attractant  near  obstructions,  and 
the  complicated  rapidly  changing  pattern  of  the  airflow. 
The  need  for  methods  to  treat  such  problems  is  demonstrated 
by  the  growing  number  of  sex  pheromone  trapping  experiments 
involving  postharvest  pests  (Sower  et  al.,  1975;  Barak  and 
Burkholder,  1976;  Read  and  Haines,  1976;  Von  Reichmuth  et  al . , 
1976,  1978;  Shapas ,  1977;  and  Vick  et  al.,  1979). 

Although  the  model  derived  below  is  applied  primarily 
to  attraction  in  a  warehouse,  it  also  applies  to  field  and 


forest  environments,  incorporating  many  elements  of  pre- 
vious models.   Like  the  Bossert  and  Wilson  (196  3)  model  it 
can  be  used  to  calculate  an  odorant  source's  attraction 
space,  a  zone  where  the  average  odorant  concentration  is 
above  a  perceptual  or  behavioral  threshold,  under  either 
still  air  or  steady  airflow  conditions.   Like  the  Aylor 
(1976)  model  it  treats  the  effects  of  turbulence  in  detail. 
Like  the  Roelofs  (1978)  model  it  includes  the  effects  of  an 
altered-behavior  threshold,  occurring  about  3  orders  of  mag- 
nitude above  the  perceptual  threshold.   Moreover,  it  con- 
siders adsorption  processes,  boundary  positions,  gravitation, 
and  the  instantaneous  structure  of  the  odorant  plume. 

The  model  is  tested  using  the  Indian  meal  moth  (IMM) , 
Plodia  inter punctel la  (Hubner)  (LepidopteraiPyralidae) .   A 
bioassay  to  determine  parameters  included  in  the  model  is 
presented  in  Chapter  I.   The  model  is  derived  and  applied  to 
several  problems  in  Chapter  II. 


CHAPTER  ONE 
ANEMOTACTIC  RESPONSE  THRESHOLD  OF  THE  INDIAN  MEAL  MOTH, 
Plodia  interparietal  la  (HUBNER)  (LEPIDOPTERA: 
PYRALIDAE) ,  TO  ITS  SEX  PHEROMONE 

An  insect's  threshold  of  response  to  sex  pheromone  is 
an  important  parameter  of  many  uses  in  olfactory  physiology 
(Kaissling,  1971)  and  applied  entomology  (see  Chapter  II) . 
The  IMM  is  a  widespread  postharvest  pest  whose  pheromonal 
physiology  has  been  studied  extensively,  but  whose  behavior- 
al response  to  pheromone  has  not  been  quantified  in  exact 
units.   IMM  populations  have  been  monitored  using  sex  phero- 
mone traps  (Vick  et  al . ,  1979);  thus,  there  is  some  prac- 
tical interest  in  the  incorporation  of  an  exactly  quanti- 
fied response  threshold  for  IMM  into  a  trapping  model.   Of 
all  the  measurable  thresholds,  the  one  most  appropriate  for 
use  in  a  trapping  model  is  the  threshold  for  upwind  anemo- 
taxis,  because  an  insect  stimulated  to  fly  has  a  greater 
likelihood  of  being  captured  by  a  trap  than  an  insect  sti- 
mulated merely  to  respond  orthokinetically  (by  antennal 
vibration  or  wing  flutter) .   To  compare  the  pheromonal  re- 
sponse of  the  IMM  with  that  of  other  insects  and  to  determine 
a  representative  threshold  for  use  in  postharvest  pest  trap- 
ping models,  the  upwind  anemotactic  response  of  the  male  IMM 
was  measured  at  23°C  and  34°C. 


Methods  and  Materials 
Insects 

Male  pupae  were  transferred  at  3-day  intervals  from  a 
laboratory  colony  (Silhacek  and  Miller,  1972)  to  an  envi- 
ronmental chamber  at  27°C  and  60%  RH  on  a  16:8  L:D  cycle 
(20-W  GE  #F201-T2-CW  light  source  in  photophase,  <0.05  lux 
in  scotophase) .   During  the  photophase,  groups  of  50-75 
newly-emerged  adults  were  placed  into  plastic  boxes  (20  x 
10  x  10  cm)  with  screen  lids,  where  they  remained  without 
food  or  water  until  testing  began  3-5  days  later. 

Pheromone 

The  IMM  pheromone,  (Z^E) -9 , 12-tetradecadien-l-ol  ace- 
tate (ZETA) ,  purchased  from  Storey  Chemical  Co.,  Willoughby, 
OH,  was  purified  twice  before  use  by  elution  with  benzene 
through  a  column  of  25%  AgNO.  on  silicic  acid.   Gas  and 
thin  layer  chromatographic  analyses  indicated  that  the 
purified  ZETA  was  at  least  99%  pure.   Between  tests  the 
ZETA  was  stored  at  -25°C. 

The  ZETA  was  dispensed  from  a  glass  tube  made  by  redu- 
cing the  unground  ends  of  a  14  cm  long  by  2  cm  ID  assembly 
of  24/40  ground-glass  joints  to  0.5  cm  ID  (Mayer,  1973). 
About  5  min  before  a  test  the  inside  of  the  assembly  was 
coated  evenly  with  a  0.5  ml  aliquot  of  ZETA  in  diethyl  ether. 
The  levels  of  ZETA  in  the  dispenser  ranged  from  0  (control) 
to  10   ng.   After  the  ether  evaporated,  the  assembled  tube 
was  stoppered  until  the  test  began. 


Olfactometers 

The  bicassay  was  done  in  3  (0.3  x  0.3  x  3.5  m)  olfac- 
tometer tunnels  described  previously  by  Mayer  (1973) .   Each 
tunnel  had  4  clcsable  sections  designated  S1-S4  in  upwind 
to  downwind  order  in  Fig.  1.   SI,  S2,  and  S4  were  56  cm  long, 
and  S3  was  140  cm  long.   The  outlet  of  a  pheromone  dispenser 
was  placed  near  the  center  of  S2  at  point  PD.   Moths  were 
released  near  the  center  of  S4  at  point  IR,  3  m  downwind 

from  PD.   The  tunnels  were  supplied  with  filtered  air  (55- 

4    3 
60%  RH,  3.6  x  10   cm  /sec  flow  rate,  35-50  cm/sec  velocity) 

in  one  of  2  temperature  ranges:   22-24  or  33-35°C.   Light 
was  provided  by  8,  rheostat-dimmed,  60  W  tungsten  bulbs 
placed  separately  inside  diffuser  boxes  spaced  uniformly 
around  the  perimeter  of  the  room  at  the  ceiling.   The  aver- 
age light  intensity  inside  the  tunnels  was  monitored  at  1 
lux  =  1.5  mw/m  . 

The  airflow  pattern  in  the  tunnels  was  depicted  by 
smoke  plumes.   The  plumes  were  produced  by  either  passing 
humidified  air  through  a  dispenser  tube  at  PD  containing  a 
cotton  swab  soaked  in  a  1:1  mixture  of  TiCl^  in  CCl,,  or  by 
substituting  a  smoke  generator  (TEM  Eng.,  Ltd.,  Crawley, 
England)  for  the  tube.   Both  methods  generated  plumes  that 
quickly  dispersed  into  less  and  less  distinct  filaments. 
Within  2  m  downwind  from  PD,  there  was  no  more  than  a  4:1 
variation  in  the  observed  smoke  density  across  the  tunnel's 
cross-section.   Because  pheromone  molecules  disperse  at 
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least  as  fast  as  smoke  particles  (Miller  and  Roelofs,  1978) , 
these  observations  indicate  that  the  average  pheromone 
concentration  at  IR  is  reasonably  approximated  by  the  ratio 
of  the  dispenser's  emission  rate  to  the  tunnel's  airflow. 

Bioassay  Procedure 

After  preliminary  experimentation  the  following  proce- 
dure was  standardized:   (1)  During  the  first  2-5  h  of  the 
scotophase,  4  0-50  male  IMM  were  loaded  into  a  screen  cage 
(17.5  x  20.3  x  25.4  cm)  and  placed  at  point  IR  in  a  tunnel. 
(2)  A  period  of  0.5-3  h  ensued  to  allow  acclimation.   (3) 
A  dosed,  stoppered  dispenser  tube  was  placed  at  PD.   The 
stoppers  were  removed  and  filtered  air  (22  C,  60  ml/min) 
was  passed  through  the  tube.   (4)  Within  10  sec  the  holding 
cage  was  opened  allowing  the  IMM  freedom  to  move  within  the 
tunnel.   After  an  additional  60  sec  the  sections  were  closed 
and  the  number  of  IMM  in  each  section  was  counted.   The 
upwind  anemotactic  response  was  estimated  by  the  fractional 
response,  Fr,  the  number  in  S2  and  S3  divided  by  the  total 
number  in  S1-S4.   (5)  The  dispenser  tubes  and  the  holding 
cages  were  washed  and  baked  for  12  h  at  200  C  after  each 
use,  and  the  tunnels  were  cleaned  with  ethanol  twice  a  week 
to  minimize  the  possibility  of  pheromonal  contamination. 

The  bioassay  had  an  incomplete  block  design.   A  treat- 
ment was  one  of  20  combinations  of  2  temperatures  (23  ,  34 

C)  and  10  tube  doses  (blank  control,  0.3,  1,  10,  102,  103, 

3     4  4  R 

3  x  10  ,  10  ,  3  x  10  ,  and  10   ng) .   Tests  were  also  done 


o  4 

at  17  C  using  10   ng  doses,  but  were  discontinued  because 

of  the  negligible  response.   A  block  was  a  set  of  6  treat- 
ments tested  on  a  given  day  at  a  given  temperature.   One 
of  the  6  was  always  a  blank,  placed  randomly  within  the 
block.   At  least  3  other  treatments  were  randomized  within 
the  remainder  of  the  block  to  eliminate  day-to-day  varia- 
tion in  the  IMM  response  from  comparison  of  the  treatments. 

Response  Analysis 

The  standard  procedure  for  calculating  a  behavioral 
threshold  is  to  adjust  the  fractional  response  for  the  con- 
trol (Abbot,  1925)  and  then  to  analyze  the  regression  of 
the  probit  adjusted  response  on  dose  (Finney,  1971) .   How- 
ever, the  stringency  of  the  upwind  anemotactic  response 
criterion  kept  the  maximal  response  in  this  bioassay  well 
below  100%,  contrary  to  the  normal  distribution  of  response 
frequency  assumed  in  probit  analysis.   If  the  response  fre- 
quency were  normally  distributed  there  would  be  a  concen- 
tration above  which  the  response  would  approach  100%,  just 
as  there  would  be  a  concentration  below  which  the  response 
would  approach  0%.   The  standard  Abbot's  correction  accounts 
only  for  deviations  from  the  normal  distribution  at  the  0% 
limit.   By  contrast,  the  following  3-step  procedure  accounts 
for  deviations  at  both  the  0%  and  the  100%  limits,  and  is 
more  generally  applicable  to  quantal  bioassays. 
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First,  the  response  was  transformed  to  probit  coor- 
dinates using  the  normalizing  equation  (Box  et  al. ,  1978): 


Nr  =  Prob  (102(Fr:Frc))  (1.1) 

Frm-Frc 


Where:   Nr     is  the  normalized  response,  in  units  of 

probits; 
Fr     is  the  uncorrected  fractional  response; 
Frc    is  the  average  fractional  response  in  the 

control; 
Frm    is  the  maximum  average  fractional  response 

at  the  given  temperature; 
Prob   is  the  integral  operator  (Finney,  1971) . 

In  the  usual  Abbot's  correction,  Frm  =  1.   Next,  the  nor- 
malized responses  from  Eq.1.1  were  fitted  by  the  standard 
probit  analysis  to  the  equation: 

Nr  =  I  +  S  Log  (DD)  (1.2) 

Where:   I      is  the  temperature-dependent  intercept,  in 
units  of  probits; 


S      is  the  temperature-dependent  slope,  probits 

-1 

ug  ; 

DD     is  the  dispenser  dose,  pg. 
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Finally  the  Nr-DD  regression  line,  Eq.  1.2,  was  converted 
to  Abbot's -corrected  coordinates  using  the  inverse  transfor- 
mation: 


p  =  ,Frm  -  FrC)  prob  1    (Nr)  (I.3) 

1  -  Frc 


Where:   Pr     is  the  Abbot's  corrected  %  response  (Abbot, 
1925); 
Prob   is  the  inverse  probit  operator. 


To  convert  behavioral  thresholds  measured  in  units  of 

\ig   dose  to  thresholds  in  the  more  practical  units  of  mole- 

3 
cules/cm  ,  the  relationship  between  a  dispenser  tube's  dose 

and  the  concentration  in  the  tunnel  must  be  quantified. 

This  relationship  was  determined  by  a  calibration  of  the 

tube's  emission  rate. 


Calibration  of  Dispenser 

At  room  temperature  pheromone  volatizes  from  a  dispen- 
ser tube  at  a  rate  proportional  to  the  initial  dose  and  in- 
versely proportional  to  the  duration  of  emission.   To  in- 
vestigate this  relationship,  emissions  from  dispenser  tubes 
loaded  with  the  pheromone  ( Z_) -7-dodecen-l-ol  acetate  (Z7AC) 
were  collected  in  1  mm  ID  x  30  cm  long  glass,  or  1  mm  ID  x 
70  cm  long,  stainless  steel  capillaries  using  a  thermal 
gradient  GLC  capillary  collector  (Brownlee  and  Silverstein, 
1968) ,  and  quantified  by  standard  gas  chromatographic  pro- 
cedures.  The  results  indicated  that  over  doses  of  1-175  yg 
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and  emission  durations  of  15-180  sec,  the  emission  rate 
was  a  constant  proportion  of  dose,  independent  of  time 
(Mayer  and  Mankin,  in  preparation) .   The  constant  of  pro- 
portionality between  the  dose  and  the  output  of  ZETA,  Ke, 
was  measured  by  a  separate  calibration  bioassay  at  23°C, 
similar  to  that  described  2  sections  previously,  except  for 
step  3  which  was  changed  to:   (31)  a  polyethylene  cap 
emitting  ZETA  at  a  rate  of  0.13  +  0.07  SE  ng/sec,  as  mea- 
sured by  the  method  of  Vick  et  al.  (197  8) ,  was  suspended  at 
PD.   The  tube  dose  evoking  the  same  average  blank-corrected 
response  as  the  cap  was  divided  into  the  measured  emission 
rate  of  the  cap  to  estimate  Ke. 

Results 
The  fraction  of  tested  IMM  attracted  to  a  dispenser 
was  sigmoidally  proportional  to  the  logarithm  of  the  dose, 
as  shown  in  Fig.  2,  and  uniformly  higher  at  34  C  than  at 
24  C.   Each  data  point  in  Fig.  2  represents  the  mean  of  8-10 
treatment  replications,  with  the  vertical  line  indicating 
the  95%  confidence  interval.   The  slope  and  intercept  para- 
meters of  each  regression  line  are  listed  in  Table  1.   Both 
the  slopes  and  the  intercepts  of  the  2  lines  are  statisti- 
cally different  (t  =  4.9  p  <0.01,  and  t  =  57.3  p  <0.001, 
respectively) .   The  good  fit  of  the  regression  lines  to  the 
data,  indicated  by  the  low  x   values  (p  >0.5  at  2  3°C  and  p 
>0.75  at  34  C) ,  supports  the  use  of  the  modified  probit  ana- 
lysis for  interpreting  pheromonal  stimulus-response  relation- 
ships. 
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The  pheromone  concentration  in  an  olfactometer  at  the 
point  of  insect  release,  IR,  was  calculated  by  assuming 
that  the  dispenser  output  dispers<-a  uniformly  in  the  3.6  x 
104  cm3/sec  airflow  (Methods:   olfactometer  section). 
The  dispenser  output  is  the  produ-'-  of  the  dispenser  dose, 
DD,  and  the  constant  of  proportio"dlitY  between  dose  and 
output,  Ke,  which  can  be  estimate^  using  point  CP  in  Fig.  1, 
the  intersection  of  the  23°C  dose  response  regression  line 
with  the  mean  Abbot's  corrected  %  response  for  8  cap  dis- 
penser test  replications.   The  Ke  ls  the  ratio  of  the  cap's 
emission  rate  to  the  tube  dose  at.  CP:   °-13  n?  sec  /2500 


ng 


=  5.2  x  10"5  sec"1.   Thus,  the  concentration  at  IR  is 


=  5.2  x  10~5  sec_1/3.6  -  10   cm   sec  ~   DD 

-9    -3 
=  1.44  x  10    cm   DD 


The  threshold  of  upwind  anem»'taxis  for  a  male  IMM  was 
calculated  using  Table  1  and  Eq.  >-4-   Following  the  stan- 
dard procedure  of  probit  analysis,  the  threshold  was  de- 
fined as  the  concentration  at  whi< !h  the  normalized  fraction 
of  response  was  0.5  (5  in  probit  •A'its)  .   Using  the  2  values 
of  DD  for  Nr  =  5  in  Eq.  1.4  gives  thresholds  of  1.34  x  10 
and  1.65  x  104  molecules  cm"3  (2.'2  x  10~   and  2-74  x  10~ 
molar)  at  23°C  and  34°C,  respectively,  which  are  compared 
with  other  reported  thresholds  in  ^ble  2. 
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Table  1.   Analysis  of  the  dose-response  regression  lines  in 
Fig.  1.   Symbols  I,  S,  Nr,  and  DD  are  defined  in 
Eqs.  1.1-2. 


Temperature 

(°C) 

Parameter 

23 

34 

I  (probits) 

5.38 

6.81 

I  Standard  error 

0.115 

0.070 

S  (probits/yg) 

0.405 

0.341 

S  Standard  error 

0.048 

0.030 

DD  for  Nr  =  5  (ug) 

3.88 

0.0048 

Upper  limit  DD  (95%) 

13.5 

0.0079 

Lower  limit  DD  (95%) 

0.054 

0.0028 

Degrees  of  freedom 

3 

4 

2 
X 

1.66  P>0.5 

1.77  P>0.75 
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Table  2.   Pheromonal  behavioral  thresholds  of  some  insects. 
Other  thresholds  are  reported  in  Kaissling  (1971) . 
The  behavioral  response  criterion  was  upwind 
anemotaxis  for  the  IMM  and  orthokinesis  for  the 
other  insects. 


Insect: 


Pheromone 


Threshold 
concentration 


(10' 


mol] 

3 

cm 


References 


Bombyx  mori: 

o 
17  C  (E,Zj-10,12-hexadecadien-l-ol    20 

21°C  10 


Kaissling  & 
Priesner  1970 


Trichoplusia  ni: 
24°C      Z7AC 

Plodia  interpunctella : 


34  C 


23  C 


ZETA 


ZETA 


16.5 
1340 


Sower  et  al , 
1971 


Chapter  I 


Trogoderma  glabrum: 

o 
27  C   (-)-14-methyl-  (Z_)  -8-hexadecenal 

2300 


Shapas  197' 


Lyman tria  dispar : 

25  C  (Z^ -7,8-epoxy-2-methyl  octadecane 

115 


Aylor  et  al, 
1976 


Discussion 
The  results  of  "the  rioassay  demonstrated  that  the  upwind 
anemotaxis  threshold  of  "the  male  IMM  for  its  sex  pheromone 
is  similar  to  the  thresholds  of  other  insects  for  their  re- 
spective sex  pheromones .   Inspection  of  Table  2  reveals  that 
the  IMM  34  C  upwind  anettiotactic  threshold  is  similar  to  the 
Bombyx  mori  (Linna^us^'.gnc*  Trichoplusia  ni  (Hubner)  ortho- 
kinetic  thresholds,  while, the  23°C  threshold  is  closer  to 
the  Trogoderma  glabrum  (Herbst)  and  Lymantria  dispar  (F.) 

orthokinetic  thresholds \      By  contrast,  human  olfactory 

8       11  3 

thresholds  generally  range  from  10   to  10   molecules/cm  , 

while  the  theoretical  limit  of  olfactory  perception  is 
about  200  molecules/c^  ".(Kaissling,  1971)  .   Because  the   4 
lowest  thresholds  in  Table  2  approach  this  limit,  it  can  be 
assumed  that  these  behavioral  thresholds  approximate  the 
corresponding  perceptual  thresholds.   In  addition,  all  of 
the  reported  ssx^pheromone  behavioral  thresholds  are  within 
4  orders  of 'l'uag'sititde^f .;#his  lower  limit.   Thus,  if  a  trap- 
ping model  is  applied  to  any  insect  whose  threshold  for 

attraction  .to  ...sex  pheronione  is  unknown,  the  threshold  can  be 

'  '■-■  3       6  3 

estimated'-. to  iije .;in .tftfea, -range  of  10   to  10  molecules/cm  . 

Dose  Dependence  of  the  Attraction  Response 

In  most  respects,  the  attraction  responses  in  Fig.  2 
follow  the  standard  sigmoidal  relationship  obtained  with 
many  other  insects  (Schneider  et  al.,  1967;  Kaissling  and 
Priesner,  1970;  Sower  et  al.,  1971;  Mayer,  1973;  Shapas,  1977) 
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In  particular,  the  change  from  0  to  100%  normalized  response 
occurs  within  about  4  concentration  decades.   However,  de- 
partures from  the  sigmoidal  relationship  appear  to  occur  at 
the  highest  tested  doses.   At  34°C  the  response  to  10  yg  is 
lower  than  the  response  to  1  yg,  and  at  23   the  response  to 
100  yg  is  lower  than  the  response  to  30  yg.   Neither  of  these 
decreases  are  statistically  significant,  but  they  are  sys- 
tematic.  In  addition,  similar  decreases  in  response  with 
increasing  dose  have  been  reported  in  olfactometer  studies 
of  other  insects  (Mayer,  1973;  Fuyama,  1976;  Hawkins,  1978), 
and  decreases  in  trap  catch  with  increasing  dose  have  been 
reported  in  field  trapping  studies  (Wolf  et  al.,  1967; 
Shorey  et  al.,  1967;  Gaston  et  al.,  1971;  Vick  et  al . , 
1979) .   This  suggests  that  the  observed  decreases  are  not 
statistical  irregularities. 

One  hypothesis  for  the  response  decrement  is  that  the 
pheromone  could  be  contaminated  with  a  slight  amount  of  the 
homolog  (Z, Z) -9 , 12-tetradecadien-i-ol  acetate,  which  is 
known  to  act  as  an  inhibitor  of  attraction  behavior  (Vick 
and  Sower,  1973) .   Such  contamination  is  unlikely,  however, 
because  the  purified  pheromone  was  found  to  be  at  least  99% 
pure  by  gas  and  thin  layer  chromatography.   More  likely,  the 
observed  decrease  is  the  effect  of  an  altered-behavior  thres- 
hold (Roelofs,  1978).   The  relationship  between  the  upwind 
anemotactic  threshold  and  the  concentration  at  which  the 
response  begins  to  depart  from  the  sigmoid  supports  this 
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hypothesis.   According  to  Roelofs  (1978),  the  altered-be- 
havior  threshold  of  an  insect  is  typically  about  3  orders  of 
magnitude  above  the  orthokinetic  threshold,  which  is  somewhat 
lower  than  the  upwind  anemotactic  threshold.   The  observed 
departure  occurs  near  10~5  ng/cm3 ,  2-3  orders  of  magnitude 
above  the  upwind  anemotactic  threshold  as  predicted. 

Temperature  Dependence  of  the  Attraction  Response 

The  effect  of  temperature  on  the  IMM  upwind  anemotac- 
tic behavioral  threshold  is  similar  to  that  reported  pre- 
viously for  the  orthokinetic  threshold  of  B.  mori  (see 
Table  2).   The  B.  mori  threshold  shifted  0.13  log  units/°C 
while  the  IMM  threshold  shifted  0.19  log  units/°C.   A  re- 
lated effect  on  the  percentage  mating  of  T.  ni  was  reported 
by  Shorey  (1966) .   The  most  likely  hypothesis  for  the  ob- 
served increase  in  threshold  with  decreasing  temperature  is 
that  the  lower  temperatures  tend  to  inhibit  flight  in  re- 
sponse to  pheromone.   Indeed,  in  several  preliminary  assays 
done  at  17  C  the  IMM  exhibited  no  response  to  10  yg  doses 
that  were  highly  stimulatory  at  23°C.   This  effect  may  be 
related  to  low  temperature  effects  on  flying  activity  re- 
ported in  a  variety  of  other  insects  (Taylor,  1963;  Bursell, 
1964;  Sanders  et  al.,  1978;  and  references  therein). 

The  findings  regarding  temperature  effects  suggest 
that  sex  pheromone  traps  may  capture  more  IMM  at  warm  tem- 
peratures than  cool  temperatures.   Similar  effects  of  tem- 
perature on  the  pheromone  trapping  of  other  insects  have 
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been  reported  by  Showers  et  al.  (1974),  Marks  (1977),  and 
Coster  et  al.  (1978) .   The  effects  of  temperature  on  trap- 
ping are  more  complicated  than  the  effects  on  attraction 
behavior,  however,  because  the  temperature  also  affects  the 
rate  of  pheromonal  emission  from  the  trap  and  temperature 
stratification  can  affect  the  pheromonal  dispersal  pattern 
(see  Chapter  II).   In  addition,  the  temperature  variation 
may  have  a  greater  influence  on  the  behavior  of  some 
insects  than  the  average  temperature  (Showers  et  al.,  1974). 

In  summary,  the  upwind  anemotactic  threshold  of  the 
male  IMM  to  its  sex  pheromone  is  similar  to  the  sex  phero- 
mone  thresholds  reported  for  other  insects.   The  threshold 
depends  on  temperature,  mostly  because  the  upwind  anemotac- 
tic response  to  pheromone  is  inhibited  at  low  temperatures. 
About  2  orders  of  magnitude  above  the  behavioral  threshold 
there  may  be  an  altered-behavior  threshold  that  causes  a 
modification  of  the  attraction  response. 


CHAPTER  TWO 
AN  INSECT  ATTRACTANCE  MODEL 


The  derivation  of  the  model  proceeds  in  stages  of  in- 
creasing complexity.   First,  the  general  problem  of  attrac- 
tion to  an  odorant  source  is  simplified  and  expressed  mathe- 
matically in  terms  of  quantifiable  parameters.   Then  the 
scope  and  complexity  of  the  problem  are  increased  by  examin- 
ing the  qualitative  effects  of  several  unquantified  factors. 
Considerable  simplification  of  the  attraction  problem  is 
achieved  by  adopting  the  following  assumptions:   (1)  An 
attractant  chemical  has  a  threshold  concentration  below 
which  its  probability  of  stimulating  an  insect  is  nil.   (2) 
An  attractant  has  a  corresponding  altered-behavior  threshold 
concentration  above  which  the  probability  of  the  insect  find- 
ing a  distant  source  decreases.   (3)  The  attraction  and 
altered-behavior  thresholds  depend  upon  the  insect's  species 
and  age,  the  temperature,  and  the  chemical  nature  of  the 
attractant,  all  of  which  are  fixed  in  a  given  application  of 
the  model.   (4)  The  emission  rate  of  the  attractant  source 
is  constant.   None  of  these  assumptions  are  strictly  valid 
under  the  usual  modeling  applications  in  which  insect  age, 
ambient  temperature,  and  source  emission  rate  vary,  but  they 
are  nevertheless  useful  heuristically . 
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Because  attraction  and  altered-behavior  thresholds  can 
be  measured  by  bioassay,  the  simplified  problem  obtained 
using  assumptions  1-4  is  essentially  solved  once  the  at- 
tractant  distribution  is  determined.   But  even  this  is  a 
formidable  task.   The  plume  emitted  from  an  attractant  source 
disperses  in  a  complicated  pattern  that  depends  on  the 
characteristics  of  the  airflow  (Skelland,  1974) .   Three  dif- 
ferent cases  will  be  considered  in  the  derivation:   dis- 
persal in  still  air  by  molecular  diffusion,  dispersal  in 
turbulent  air  currents  of  zero  average  speed  and  direction, 
and  dispersal  in  turbulent  or  laminar  air  currents  of  con- 
stant average  speed  and  direction. 

Case  I:   Molecular  Diffusion 
In  still  air,  an  attractant  plume  disperses  by  molecu- 
lar diffusion,  which  is  described  by  the  mass-balance  equa- 
tion (Veigele  and  Head,  197  8)  : 


q  =  (-ft-  _Dv2)  C,  (II. 1) 


Where:   q      is  the  rate  of  emission  per  unit  volume,  with 

units  of  g  sec-   cm"  ; 

2     —1 

D      is  the  diffusion  coefficient,  cm   sec   ; 

C      is  the  attractant  concentration,  g  cm   ; 


_3_ 


,2 


represents  differentiation  with  respect  to  time; 
represents  the  Laplacian  differential  operator 
(Protter  and  Morrey,  1966;  p.  567). 


24 


Theoretically  Eq.  II. 1  has  many  solutions,  depending  on  the 
initial  concentration  distribution  and  the  first  partial 
derivative,  3C/3r,  at  all  boundary  surfaces.   For  heuristic 
purposes  it  is  convenient  to  assume  that  the  source  is  sur- 
rounded by  a  single  boundary  surface  and  the  initial  (t  =  0) 
concentration  is  zero  everywhere  inside  the  boundary.   The 
value  of  3C/3r  at  this  boundary  is  (Chamberlain,  1953; 
Judeikis  and  Stewart,  1976;  Draxler  and  Elliot,  1977): 

V 
3C  _  d 
37  ~  D~  C'  (II. 2) 

Where:   V^     is  an  empirical  parameter,  the  deposition 
velocity,  cm  sec-1. 

The  solution  to  Eqs.  II. 1-2  for  dispersal  inside  a  spheri- 
cal boundary  is  (Carslaw  and  Jaeger,  1967,  p.  367,  and  Ap- 
pendix) : 

2     2   2 
(ah-1)   +  a   9 

Q  n 

u    2TrarD  L         2  2 

n=l  a  e    +  ah (ah-1) 
n 

[sin(re  )/9  ]  [l-exp(-D62t)  ],       (II. 3) 
n   n  n 

Where:   Q      is  the  emission  rate  of  the  source,  g  sec"  ; 

a      is  the  radial  distance  of  the  boundary  sphere, 
from  the  source  cm; 

r      is  the  radial  distance  of  the  measurement  po- 
sition from  the  source,  cm; 
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t    is  the  duration  of  emission,  sec; 

h    is  the  ratio  V  /d,    cm   ; 

8n   is  the  nth  positive  root  of  the  equation 

a  cot  (a9)  +  ah  =  1.  (II. 4) 

Equations  II. 3-4  are  not  restricted  particularly  to 
attraction  problems  with  spherical  boundaries  because  most 
attractants  have  values  of  D  and  V,  that  fall  within  fairly 
narrow  ranges,  limiting  the  extent  of  boundary  effects.   The 

molecular  diffusion  coefficient  of  a  sex  pheromone  is  about 

2 
0.05  cm  /sec  (Wilson  et  al.,  1969;  Hirooka  and  Suwanai,  1976), 

and  most  other  attractants  have  molecular  diffusion  coeffi- 

2 
cients  near  0.03-0.07  cm  /sec  (Monchiek  and  Mason,  1961; 

Lugg,  1968) .   However,  unless  conditions  are  strictly  con- 
trolled, air  currents  usually  occur  that  effectively  in- 

2 
crease  D  to  0.1-0.5  cm  /sec  (Bossert  and  Wilson,  1963).   The 

magnitude  of  the  deposition  velocity  depends  primarily  upon 
the  forces  causing  adsorption  of  the  attractant  vapor  to 
the  substrate.   Interfaces  composed  of  different  vapors  and 
different  substrates  have  similar  deposition  velocities  inas- 
much as  the  vapor-substrate  interactive  forces  are  similar 
(Judeikis  and  Stewart,  1976) .   Measurements  of  deposition 
velocity  typically  vary  over  the  range  0.1  to  10  cm/sec  for 
Iodine-plastic,  SO  -concrete,  and  pheromone-vegetation  in- 
terfaces (Chamberlain,  1953,  1966;  Judeikis  and  Stewart, 
1976;  Nakamura  and  Kawasaki,  1977).   The  following  analysis 
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Of  Eqs.  II. 3-4  shows  that,  when  D  and  V"d  fall  within  the 

ranges  given  above,  the  effect  of  a  boundary  on  the  attrac- 

tant  distribution  is  usually  small. 

The  spatial  variation  of  the  relative  concentration, 

Cr  =  C/Q,  is  depicted  in  Figs.  3-6  at  several  different 

values  of  a,  D,  V  ,  and  t.   The  curves  were  calculated 
d 

from  the  first  700  terms  of  the  series  in  Eq.  II. 3,  which 
converged  after  about  200  terms  except  at  small  r  and  t. 
The  curve  for  a.  =  "~°  corresponds  to  the  boundless  case,  pre- 
viously discussed  by  Bossert  and  Wilson  (1963) ,  whose  solu- 
tion reduces  to: 


Cr  =   o-T^r-   erfc  (r/  (4Dt) 1/2)  ,  (II.  5) 


Where:   erfc    is  the  complimentary  error  function  (Carslaw 
and  Jaeger,  1967) . 

Inspection  of  Figs.  3-4  shows  that  the  smaller  the 

boundary  radius  the  smaller  the  variation  of  C  with  time, 

and  as  t  increases  the  difference  between  C   in  a  bounded 

r 

and  a  boundless  environment  decreases.   The  2  curves  for 
a  =  °°  differ  considerably  from  each  other,  but  the  curve 
for  a  =  150  at  t  =  60  sec  is  quite  similar  to  the  curve  for 
a  =  100  at  t  =  8.6  x  105  sec.   After  about   103  sec,  there 
is  very  little  difference  between  the  curves  a  =  °= ,  a  =  1000, 
and  a  =  100,  short  distances  away  from  the  respective  boun- 
daries.  The  curves  for  103  sec  are  not  plotted  because  they 


Figure  3.   Spatial  variation  of  the  relative  concentration, 
Cr  =  C/Q,  at  t  =  60  sec.   The  regression  lines 
are  calculated  from  Eqs.  II. 3-4  with  D  =  1  cm2/ 
sec,  and  Vd  =  1  cm/sec.   The  radius  of  the  boun- 
dary sphere,  a,  has  units  of  cm. 
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Figure  4.   Spatial  variation  of  the  relative  concentration, 
Cr  =  C/Q  at  t  =  8.6  x  105  sec.   The  regression 
lines  are  calculated  from  Eqs .  II. 3-4  with  D  =  1 
cm  /sec,  and  Vd  =  1  cm/sec.   The  radius  of  the 
boundary  sphere,  a,  has  units  of  cm. 
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Figure  5.   Variation  of  the  relative  concentration,  C  = 

C/Q,  with  respect  to  the  diffusion  coefficient  D. 
The  regression  lines  are  calculated  from  Eqs . 
II. 3-4  with  t  =  8.6  x  105  sec,  Vd  =  1  cm/sec, 
and  a  =  1000  cm.   D  has  units  of  cm2/sec. 
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Figure  6.   Variation  of  the  relative  concentration,  Cr  - 
C/Q,  with  respect  to  the  deposition  velocity, 
Vj.   The  regression  lines  are_calculated  from 
Eqs.  II. 3-4  with  t  =  8.6  x  105  sec,  D  =  1  cm2/ 
sec,  and  a  =  1000  cm.   V-,  has  units  of  cm/sec. 
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are  quite  similar  to  the  curves  in  Fig.  4.   Indeed,  the 
effect  of  the  boundary  can  be  disregarded  so  long  as  r/a< 

0.9,  unless  D  is  considerably  smaller  than  1  cm  /sec. 

_2 
This  is  shown  in  Figs.  5-6,  where  D  ranges  from  10    to 

2    2  -5       2 

10   cm  /sec,  and  Vd  from  10    to  10   cm/sec,  at  t  =  8 . 6  x 

10   sec  and  a  =  1000  cm.   Even  at  the  hypothetical  lower 

2 
limit  of  the  molecular  diffusion  coefficient,  0 . 1  cm  /sec, 

the  influence  of  the  boundary  is  negligible  until  r/a>0.75 

Under  these  conditions,  the  exact  geometry  of  a  boundary 

is  unimportant,  and  the  solution  to  Eqs.  II. 1-2  under 

any  boundary  geometry  is  similar  to  Eqs.  II. 3-4.   This 

justifies  the  general  use  of  Eqs.  II. 3-4  in  a  model  for 

attraction  in  still  air,  provided  that  the  results  are 

interpreted  with  caution  and  when  the  emission  duration 

is  less  than  about  an  hour  and/or  r/a>0.9. 

Case  II:   Dispersal  in  Airflow  of  Zero  Average  Velocity 
The  next  step  of  the  derivation  is  to  consider  the 
effect  of  random  air  currents  on  the  diffusion  process. 
In  most  applications  of  a  model  the  air  has  some  movement. 
Often,  random  fluctuations  of  speed  and  direction  occur, 
caused  by  superimposed  whirls  or  eddies  of  various  sizes. 
Eddies  with  a  diameter  greater  than  about  100  cm  or  less 
than  about  5  cm  have  little  effect  on  an  attractant  plume 
but  eddies  of  intermediate  diameter  cause  the  plume  to 
disperse  rapidly  (Aylor,  1976) . 
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This  kind  of  dispersal,  called  turbulent  or  eddy  dif- 
fusion, is  analogous  to  molecular  diffusion.   If  the  air- 
flow has  a  near-zero  average  velocity,  eddy  diffusion  is 
described  via  Eq.  II. 1,  replacing  the  molecular  diffusion 
coefficient  with  the  sum  of  the  molecular  and  eddy  diffu- 
sion coefficients.   The  eddy  diffusion  coefficient  of  a 

o 

vapor  or  an  aerosol  is  0.1-10  cm  /sec  in  a  calm  environment 

(Sutton,  1953;  Pasquill,  1961;  Bossert  and  Wilson,  1963; 
Allen,  1975).   Just  as  in  molecular  diffusion,  Eqs .  II. 3-4 
can  be  used  as  a  general  solution  to  the  turbulent  diffu- 
sion problem,  provided  that  the  results  are  interpreted 
with  caution  when  t<3600  sec  and/or  r/a>0.9. 

Case  III:   Dispersal  in  Airflow  of 
Constant  Average  Velocity 

The  derivation  of  an  attraction  model  is  now  completed 
by  considering  the  effect  of  relatively  stable  convection 
on  molecular  and  turbulent  diffusion  processes.   An  exact 
mathematical  description  of  convective-dif fusive  dispersal 
does  not  exist,  but  many  time-averaged  statistical  descrip- 
tions have  been  published  (e.g.,  Sutton,  1953;  Pasquill, 
1961;  Lamb  et  al . ,  1975).   The  Sutton  equation  is: 


C  =  —iQ expt-X1-75  (iL  +  ii,)],  (II'6) 

1.75  K  z    k 

K  K  vX  y      y 

y  z  y 
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where:   C      is  the  attractant  concentration,  with  units 
of  g  sec   ; 

Q      is  the  attractant' s  emission  rate,  g  sec"  ; 

v      is  the  mean  air  velocity,  cm  sec   ; 

X      is  the  distance  from  the  source  along  the 
airflow  axis,  cm; 

Y      is  the  distance  from  the  source  along  cross- 
wind  horizontal  axis,  cm; 

Z      is  the  distance  from  the  source  along  the 

vertical  axis,  cm; 

K  ,K   are  empirical  constants  whose  values  have  been 
Y   z 

tabulated  for  different  turbulence  levels, 
cm0'125  (Gifford,  1960). 

Equation  II. 6  has  been  incorporated  into  several  of  the  mo- 
dels cited  in  the  introduction. 

However,  there  are  several  cautions  to  using  Eq.  II. 6. 
When  barriers  or  thermal  stratification  interfere  with  free 
air  movement,  or  the  airflow  becomes  highly  unstable,  tab- 
ulated values  of  K   and  K  are  not  very  reliable  (Sutton, 
1953).   Under  such  conditions,  the  dispersal  pattern  can  be 
observed  directly  by  using  smoke  plumes,  and  the  tabulated 

values  of  K  and  K  can  be  replaced  with  values  determined 
y      z 

by  the  extent  of  the  plumes  (Pasquill,  1961) .   Even  allowing 
for  these  effects,  Eq.  II. 6  provides  no  information  about 
the  instantaneous  concentration  distribution,  which  is  a 
more  important  determinant  of  the  insect's  searching  behavior 
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than  the  average  distribution  (Aylor,  1976) .   Qualitative 
effects  of  differences  in  the  instantaneous  distribution  can 
nevertheless  be  incorporated  into  the  model  as  follows. 

An  attractant  plume  typically  consists  of  3  distinct 
regions  that  differ  in  shape  (Aylor,  1976;  J.  Kittredge, 
personal  communication).   Region  I  extends  0.5-10  m  down- 
wind of  the  source.   The  plume  in  this  region  is  a  single 
filament  whose  width  is  smaller  than  the  widths  of  the 
smallest  eddies;  thus,  it  disperses  solely  by  the  slew  pro- 
cess of  molecular  diffusion.   It  splits  into  several  fila- 
ments in  Region  II,  which  extends  from  the  edge  of  Region  I 
to  about  50-70  m  downwind,  depending  on  the  turbulence  and 
the  air  velocity  (Hinze,  1975) .   In  addition,  it  begins  to 
expand  rapidly  because  its  width  is  now  comparable  to  the 
widths  of  the  smaller  eddies,  allowing  dispersal  by  tur- 
bulent diffusion.   This  expansion  continues  until  the  limits 
of  the  plume  become  indistinct  in  the  final  region,  which 
extends  downwind  indefinitely  from  the  edge  of  region  II. 
The  plume  has  greater  definition  in  region  I  than  in  region 
II,  and  the  least  definition  in  region  III. 

The  regional  differences  in  plume  shape  can  be  used  to 
improve  the  model's  estimates  of  the  attraction  space.   Ac- 
cording to  current  theories  of  flight  orientation  behavior 
(Farkas  and  Shorey,  1974;  Kennedy,  1974;  and  references 
therein) ,  an  insect  finds  it  easier  to  steer  toward  a  source 
when  the  plume  structure  is  well-defined,   A  model  incorpor- 
ating this  hypothesis  predicts  that  an  insect  is  more  likely 
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to  reach  the  source  if  it  enters  region  I  than  if  it  enters 
the  other  2  regions.   A  corollary  is  that  the  smaller  the 
dimensions  of  the  source,  the  greater  the  extent  of  region 
I,  and  the  greater  the  probability  of  an  insect  finding 
the  source . 

The  choice  whether  to  use  Eq.  II. 6  or  Eqs .  II. 3-4  in 
a  particular  application  of  the  model  depends  upon  the  tur- 
bulence level,  the  obstructing  boundaries,  and  the  shape  of 
the  plume.   If  the  shape  of  the  plume  is  unknown,  it  can  be 
determined  by  smoke  plume  observations  (see  Chapter  I) . 
Generally  Case  I  dispersal,  characterized  by  the  absence  of 
turbulence,  occurs  only  in  a  highly  controlled  environment. 
If  a  plume  in  turbulent  air  is  nearly  spherical,  Case  II  ap- 
plies, and  if  it  is  ellipsoidal,  Case  III  applies.   The  move- 
ment of  a  plume  in  sheltered  areas,  cul-de-sacs,  or  a 
closed,  empty  warehouse  with  an  isothermal  temperature  dis- 
tribution usually  fits  a  Case  II  dispersal  pattern,  while 
Case  III  dispersal  is  more  likely  to  occur  in  open  venti- 
lated corridors  or  highly  stratified  air. 

Before  proceeding  with  applications  of  the  model  it  is 
appropriate  to  consider  briefly  the  validity  of  Assumptions 
1-4  and  the  precision  of  the  dispersal  equations.  The  sim- 
plified problem  obtained  by  adopting  these  assumptions  ne- 
glects such  behavioral  factors  as  visual  attraction  (Shorey 
and  Gaston,  1965;  Hienton,  1974),  anemotaxis  (Kennedy,  1974), 
habituation  (Thompson  and  Spencer,  1966;  Traynier,  1968; 


40 


Sower  et  al . ,  1973;  Marks,  1978)  and  changes  in  the  in- 
tensity of  a  response  at  different  attractant  concentrations 
(Cain  and  Engen,  1969;  Mayer,  1973;  Bartell  and  Lawrence, 
1977) .   None  of  these  factors  are  quantified  very  precisely; 
consequently,  their  inclusion  in  the  model  would  not  neces- 
sarily improve  its  precision  at  this  time.   The  model  also 
neglects  such  physical  factors  as  pheromonal  lability  (Lund- 
berg,  1961;  Sower  et  al . ,  1975)  differences  in  the  ratios 
of  attractant  components  (Roelofs,  1978),  and  changes  in 
the  emission  rate  of  the  source  at  different  temperatures 
and  airflows.   These  factors,  as  well  as  the  statistical 
nature  of  atmospheric  mass-transfer  processes,  limit  the 
precision  of  the  dispersal  equations  used  in  the  model. 
Because  of  these  limitations,  the  model  must  be  used  in 
conjunction  with,  rather  than  instead  of  experimental 
studies.   Notwithstanding,  the  model  provides  considerable 
insight  into  the  attraction  process,  as  shown  by  the  follow- 
ing applications. 

Results 
The  applications  presented  here  are  concerned  primari- 
ly with  Case  II  dispersal  because  Case  I  is  of  little  prac- 
tical interest  and  Case  III  has  been  treated  in  detail  by 
the  models  cited  earlier.   Under  Case  I  dispersal,  the  model 
treats  the  hypothesis  (Traynier,  1968;  Perez  and  Hensley, 
1973)  that  pheromone-laden  air  tends  to  sink.   Under  Case  II, 
attraction  spaces  and  altered-behavior  spaces  are  calculated 
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for  calling  insects  and  sex  pheromone  traps.   The  model 
also  deals  with  attraction  to  competing  sources.   Under 
Case  III,  the  model  entertains  the  possibility  that  insects 
can  be  attracted  from  outside  a  warehouse. 

Case  I  Dispersal:   Effect  of  Gravity  on  Pheromone  in  Still 
Air 

Equations  II. 1,  3,  and  4  do  not  consider  gravitation; 
thus,  by  default  they  predict  that  the  dispersal  of  phero- 
mone in  still  air  is  unaffected  by  the  high  molecular 
weight  of  pheromone  molecules  relative  to  the  weight  of  air 
molecules.   It  will  now  be  shown  that  this  prediction  re- 
mains unchanged  after  gravitation  is  incorporated  into  the 
model. 

The  effect  of  gravity  is  determined  by  Archimedes' 
principle  (Sears,  1958,  p.  365),  in  that  an  air-pheromone 
mixture  is  subject  to  a  gravitational  force  proportional  to 
the  difference  between  the  density  of  the  mixture  and  the 

density  of  the  air  surrounding  it.   The  standard  density  of 

.   .  -3      3 

air  is  p   =  1.3  x  10   g/cm   (Sears,  1958).   The  density  of 
a 

the  air-pheromone  mixture  is,  by  definition: 


pa   =  m  +  m  (II. 7) 

ap   _a    p 


Where:   m      is  the  mass  of  air  inside  V,  g; 
a 

-  m      is  the  mass  of  pheromone  inside  V,  g; 
P 

V      is  the  volume,  cm  . 
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The  pheromone  vapor  density,  m  /V,  reaches  a  maximum 
at  partial  pressures  approaching  the  vapor  pressure,  the 
pressure  exerted  by  the  pheromone  vapor  when  the  air  is 
saturated.   The  vapor  pressure  and  the  corresponding  sa- 
turated vapor  density  of  a  12-16  carbon  sex  pheromone  are 

-4  3 

about  1  x  10    cm-Hg  and  10  ng/cm  ,  respectively  (Hirooka 

and  Suwanai,  1976).   The  magnitude  of  m  /V  in  Eq.  II. 7  can 
be  calculated  by  combining  Dalton's  Law  of  Partial  Pres- 
sures, the  Ideal  Gas  Law,  and  the  relationship  between 
mass  and  molecular  weight  which  are  respectively, 


P  =  ^  =  Pr,  +  P  '  (II. 8) 

at    P    a 


PV  =  NRT,  (II. 9) 


and 


m  =  NM,  (11.10) 


Where:   P      is  the  pressure  in  units  of  cm  Hg  or  dyne 

-2 
cm   ; 

N      is  the  number  of  moles; 

7  -1 

R      is  the  gas  constant,  8.31  x  10   ergs  mole 

Kelvin; 

T      is  the  temperature,  deg  Kelvin; 

m      is  the  mass,  g; 

M      is  the  molecular  weight,  g/mole. 
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The  subscripts  a,  p,  and  at,  refer  to  unadulterated 
air,  pheromone,  and  atmosphere,  respectively.   The  result 
of  combining  Eqs .  II. 8-10  is 


m  /V  =  p   -  m  M  /VM  .  (11.11) 

a      a    p  a   p 


Accordingly,  Eq.  II. 7  can  be  rewritten: 


P    =  P   +  m   (1-M  /M  )/V.  (11.12) 

aP    a    p     a   p 


M  ,  the  molecular  weight  of  air,  is  about  29  g/mole  and  M  , 
a  p 

the  molecular  weight  of  the  pheromone,  is  about  225  g/mole. 

3 
Given  a  maximum  value  of  10  ng/cm   for  m  /V,  the  maximum 

density  of  the  air-pheromone  mixture  is: 


-3  -9      3 

P    =  1.3  x  10    +  8.7  x  10    g/cm  .  (11.13) 
ap 


This  shows  that  there  is  essentially  no  difference  between 
the  density  of  the  air-pheromone  mixture  and  the  unadulter- 
ated air;  so  by  Archimedes'  principle,  gravity  has  little 
effect  on  the  mixture. 

By  contrast,  a  solvent  such  as  diethyl  ether  may  sink 
in  still  air  because  it  has  a  high  vapor  pressure,  about 
442  mm-Hg.   Using  Eq.  II. 9  and  the  conversion  factor  1  cm-Hg 
=  1.33  x  10   dynes/cm, 
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-3   ,   3 
m  =  M  P  =  1.75  x  10    g/cm  ,         (11.14) 
e    e 

V     RT 


where:   M      is  the  molecular  weight  of  diethyl  ether, 
e 

74  g/mole; 

m      is   the  mass  of  ether  inside  V,  g. 
e 


Accordingly,  the  density  of  the  air-diethyl  ether  mixture 

is  P    =  2.36  x  10~3  g/cm3. 
ae 

The  net  force  on  1  ml  of  ether-saturated  air  is,  by 
Archimedes'  principle: 


4     2 
F  =  (p    -p  )  980  cm  /sec   =  1.04  dynes. 
ae   a 

(11.15) 

After  1  sec  a  ml  of  ether-saturated  air  falls  at  the 
velocity 

F 
V  =  "o x  1  ml  (1  sec)  =  440-7  cm/sec. 

ae 

(11.16) 

This  is  about  half  the  rate  of  fall  of  a  solid  object.   How- 
ever, immediately  after  the  mixture  is  emitted  it  begins  to 
disperse  into  the  air.   Within  a  short  period  of  time  the 
density  inside  V  equilibrates  with  that  of  the  surrounding 
air  ana  the  rate  of  fall  decreases  to  zero. 
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Case  II  Dispersal :   Attraction  and  Altered  Behavior 
Spaces  of  IMM 

If  the  attraction  and  altered-behavior  thresholds  of 

an  insect  are  known,  Figs.  3  and  4  provide  estimates  of  the 

attraction  and  altered-behavior  spaces  of  an  attractant 

source  in  a  Case  II  environment.   The  insect  considered 

here  is  the  widespread  postharvest  pest,  P.  intcrpunctella 

(IMM) .   A  male  IMM  has  an  upwind  anemotactic  threshold  of 

6.8  x  10    ng/cm   (10     molar)  at  34  C,  5.6  x  10    ng/cm 

at  23  C,  and  it  has  a  hypothetical  altered-behavior  thres- 

-5       3 
hold  of  about  10    ng/cm   (see  Chapter  I) .   A  calling  vir- 

-4 
gin  female  IMM  emits  pheromone  at  8  x  10    ng/sec  (Sower  and 

Fish,  1975) ,  and  typical  sex  pheromone  traps  for  capturing 
IMM  emit  0.01-0.76  ng/sec  (Vick  et  al . ,  1979).   The  rela- 
tive thresholds  for  each  emission  rate  and  temperature  are 
given  in  Table  3.   Insertion  of  the  IMM  female's  relative 
attraction  threshold  into  Fig.  3  demonstrates  that  the  at- 
traction space  of  a  female  after  one  min  of  calling  is 
rather  insensitive  to  temperature  but  somewhat  dependent 
on  boundary  position.   The  entire  volume  inside  a  boundary 
of  150  cm  radius  is  an  attraction  space.   The  attraction 
space  inside  a  boundary  of  1000  cm  radius  is  a  sphere  of 
about  250  cm  radius.   In  a  boundless  environment,  the  at- 
traction sphere  has  a  radius  of  about  4  0  cm.   Actually, 
there  is  little  practical  difference  among  these  radii.   The 
spaces  calculated  by  the  model  for  the  IMM  are  compared  with 


4  6 


Table  3.   Relative  attraction  thresholds,  C   =  C/Q,  and 

relative  altered-behavior  thresholds,  Crcj ,  for  a 
female  I MM  and  2  sex  pheromone  traps  emitting  at 
different  rates.   Units  of  C   and  Cr(j  are  sec/ 
cm3. 


Female 

Pheromone 

Pheromone 

Parameter 

IMM 

trap 

trap 

Q  (ng/sec) 

8  x  10~4 

0.01 

0.76 

C   at  23°C 

r 

-4 
7  x  10 

5.6  x  10~5 

7.4  x  10~7 

34°C 

8.5  x  10" 

-6 

6.8  x  10~7 

8.9  x  10-9 

C 
rd 

1.25  x  10" 

•2 

lO"3 

1.3  x  10~5 
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the  spaces  for  different  insects  in  Table  4.   The  difference 
between  the  first  2  estimates,  derived  from  Eq.  II. 1,  and 
the  last  2  estimates,  derived  from  Eq.  II. 6,  indicates 
primarily  the  effect  of  convection  on  the  pheromone  disper- 
sal pattern. 

Hitherto,  the  concept  of  an  altered -behavior  space  has 
been  applied  to  traps  but  not  to  calling  insects,  perhaps 
because  it  had  been  assumed  that  the  concentration  around  a 
calling  insect  does  not  rise  above  the  altered-behavior 
threshold.   According  to  the  model  however,  calling  females 
of  both  IMM  and  Trichoplusia  ni  (Hubner)  species  are  sur- 
rounded by  altered  behavior  spaces  of  small  but  definite 
extent.   Using  Table  3  and  Fig.  3,  the  altered-behavior 
threshold  of  a  female  IMM  occurs  at  about  6  cm.   A  female 
Trichoplusia  ni  (Hubner)  emits  pheromone  at  the  rate  of  0.1 
ng/sec,  and  a  male  has  an  activation  threshold  of  about 

3  x  10   ng/cm   (Sower  et  al.,  1971).   If  it  is  assumed  that 

-4      3 
the  altered  behavior  threshold  is  10    ng/cm  ,  3  orders  of 

magnitude  above  the  activation  threshold  (Roelofs,  1978), 

-3        3 
the  relative  threshold  is  C  ,  =  10    sec/cm  ,  and  from  Fiq. 

rd  3 

3,  the  altered-behavior  space  is  a  sphere  of  about  15-60 
cm  radius. 

The  attraction  and  altered-behavior  spaces  of  IMM  sex 
pheromone  traps  can  be  determined  using  Table  3  and  Fig.  4. 
If  either  trap  in  Table  3  is  the  source,  essentially  the 
entire  volume  inside  a  boundary  of  10  or  1000  cm  radius  is 
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Table  4.   Calculated  values  for  the  maximal  communication 
distance  or  attractive  range  of  an  insect,  using 
either  Eqs.  II. 3-4  or  Eq.  II. 6. 


Insect 


Range (m) 


Reference 


Plodia  interpunctella  (Hubner) 
Pogonomyrmex  badius  (Latreille) 

Trogoderma  glabrum  (Herbst) 
Hyphantria  cunea  (Drury) 

Spodoptera  litura  (F . ) 

Trichoplusia  ni  (Hubner) 


0.4-2.5   Chapter  II 

1.04      Bossert  &  Wilson 

(1963) 
0.8-10    Shapas  (1977) 
3-10      Hirooka  & 

Suwanai  (1976) 
80       Nakamura  & 

Kawasaki  (1976) 
1-100     Sower  et  al. 

(1971) 
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above  the  attraction  threshold.   At  0.76  ng/sec  the  entire 
volume  is  also  above  the  altered-behavior  threshold.   By 
contrast,  the  0.01  ng/sec  trap  has  an  altered-behavior 
sphere  of  71  cm  radius  inside  a  1000  cm  boundary  and  45  cm 
radius  inside  a  100  cm  boundary.   This  suggests  that  a  trap 

emitting  1  ng/sec  might  capture  fewer  I MM  than  an  otherwise 

-3 
identical  trap  emitting  10    ng/sec. 

Because  the  male  IMM  anemotactic  threshold  to  sex 
pheromone  is  about  2  orders  of  magnitude  higher  at  23  C 
than  at  34  C,  it  is  surprising  that  there  is  little  effect 
of  temperature  on  the  predicted  attraction  space  of  a  cal- 
ling female  or  a  pheromone  trap.   This  result  is  due  to  a 

rapid  decrease  in  C  with  respect  to  position  once  C   falls 
r  r 

-4       3 
below  about  10    sec/cm  .   The  temperature  could  have  a 

-4 
large  effect  if  the  threshold  C  were  higher  than  10    sec/ 

3 
cm  .   Such  high  threshold  values  of  C  may  occur  for  non- 

pheromonal  attractants . 

The  model  cannot  be  used  to  calculate  attraction  to 

competing  sources  unless  the  attraction  spaces  are  nonover- 

lapping  or  the  principal  searching  mechanism  is  chemotaxis. 

If  the  insect  searches  by  additional  mechanisms,  e.g. 

klinokinesis ,  chemokinesis  or  vision  (Wright,  1958;  Farkas 

and  Shorey,  19  74;  Wail  and  Perry,  197  8;  Baker  and  Carde, 

1979)  the  problem  of  attraction  to  competing  sources  becomes 

very  complicated,  particularly  if  trap-female  competition  is 

considered.   For  example,  a  female  IMM  generally  calls  from 
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walls,  ceilings,  or  other  exposed  surfaces.   Typically,  a 
male  I MM  stimulated  by  sex  pheromone  orients  visually  to 
such  surfaces,  investigating  objects  resembling  female  IMM 
(Sower  et  al.,  1975).   This  behavior  increases  the  proba- 
bility of  locating  a  female  but  decreases  the  probability 
of  capture  by  a  trap  unless  the  trap  is  highly  attractive 
visually.   Thus,  as  the  density  of  female  IMM  increases, 
a  trap  could  lose  efficiency  much  more  rapidly  than  the 
model  would  predict  from  camparisons  of  the  attraction 
spaces . 

However,  if  the  attraction  spaces  are  nonoverlapping 
or  the  principal  searching  mechanism  is  chemotaxis,  attrac- 
tion to  competing  sources  can  be  calculated  by  superposi- 
tion of  individual  solutions  of  the  model  equations.   The 
procedure  for  nonoverlapping  sources  has  been  treated  in 
detail  elsewhere  (Knipling  and  McGuire,  1966;  Nakamura  and 
Oyama,  1978) ,  so  only  attraction  under  chemotaxis  will  be 
considered  here.   Suppose  2  sources,  El  and  E2,  are  located 
500  cm  apart  in  a  warehouse  whose  length-width-height  dimen- 
sions are  about  20  m  each.   If  the  emission  rates  of  El  and 
E2  are  equal  and  the  attraction  spaces  overlap,  an  insect 
inside  the  overlapping  region  will  fly  to  the  closest  source, 
If  the  emission  rates  are  unequal,  the  probability  of  at- 
traction to  either  source  can  be  predicted  using  Fig.  4, 
which  indicates  that  a  10-fold  decrease  in  pheromone  concen- 
tration is  roughly  equivalent  to  a  10-fold  increase  in  the 
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distance  from  the  source.   Accordingly,  when  the  emission 
rate  of  El  is  10-fold  greater  than  that  of  E2 ,  an  insect 
inside  the  overlapping  space  flies  to  El  unless  it  is 
within  about  4  5  cm  of  E2. 

Case  III  Dispersal:   Attraction  of  Insects  from 
Outside  a  Warehouse  

Equations  II. 3-4  can  be  combined  with  Eq.  II. 6  to 
determine  whether  emission  from  a  trap  inside  a  warehouse 
produces  an  above-threshold  concentration  of  pheromone  out- 
side.  Suppose  that  a  trap  emitting  IMM  sex  pheromone  at  the 
rate  of  0.1  ng/sec  is  placed  near  the  middle  of  a  warehouse 
whose  length-width-height  dimensions  are  about  20  m  each. 
After  several  days  a  2  x  2  m  door  is  opened.   The  ambient 
temperature  is  27°C,  and  the  outside  airflow  is  a  constant 
50  cm/sec.   The  Ky  and  K   in  Eq.  II. 6  are  assumed  to  have 
the  values  tabulated  by  Sutton  (1953)  for  dispersal  under 
stable  conditions,  0.4  and  0 . 2  cm  *    ,  respectively.   The 
maximum  rate  of  flow  of  pheromone  through  the  door  is 
(Judeikis  and  Stewart,  1976): 


Q  -  (RT/2itM  )1/2  C(4  x  104  cm2).        (11.17) 


The  values  of  R  and  M   are  given  in  the  Results:   Case  I 
section.   The  concentration,  C,  is  found  by  inspection  of 

Fig.  4,  which  indicates  that  the  pheromone  concentration 

—  fi       3  —  5        3 

near  the  door  is  about  10   ng/cm   (C  =  10    sec/cm  ) . 
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Inserting  these  values  into  Eq.  11.17  gives  Q  =  168  ng/sec. 
This  rate  decreases  immediately  after  the  door  is  opened, 
and  continues  to  decrease  to  0.1  ng/sec.   For  estimation 
purposes  it  is  convenient  to  choose  Q  =  10  ng/sec  and  Y  = 
Z  =  0  in  Eq.  II. 6,  which  then  reduces  to 


-9      -1.25   -1.75 
C  =  1.59  x  10    g  cm      x      .     (11.18) 


-17      3 
The  behavioral  threshold  of  a  male  IMM  is  ca .  10    g/cm 

(see  Chapter  I) .   Accordingly,  the  maximum  downwind  dis- 
tance from  which  a  male  IMM  can  be  attracted  is 


X  =  (10-17  g  cm  3/1.59  x  10  9  g  cm  1*25'  '0-571 
-  4.82  x  104  cm.  (11.19) 


2 

If  the  pheromone  is  emitted  from  a  1  cm  hole  instead 

of  a  2  m  x  2  m  door,  the  emission  rate  is  Q  -  4.2  x  10 
g/sec.   At  this  rate  Eq.  II. 6  reduces  to 


,-."17      -3.,  rn         ,  -13      -1.25,  -0.571 
X  =  (10     g  cm  /6. 68  x  10     g  cm 

=  569  cm.  (11.20) 


Thus,  the  attraction  of  insects  from  outside  the  warehouse 

cannot  be  neglected  unless  the.  external  openings  are  less 

2 

than  about  1  cm   area. 
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Discussion 

Some  of  the  predictions  resulting  from  the  incorpora- 
tion of  boundary  effects,  plume  shape  effects,  altered- 
behavior  thresholds,  and  gravitation  into  the  attraction 
model  warrant  further  discussion.   These  factors  are  con- 
sidered in  order  below,  in  the  context  of  previously  pub- 
lished work. 

Inclusion  of  boundary  effects  in  the  model  results  in 
2  predictions.   First,  according  to  Figs.  3-4,  the  position 
of  the  boundary  is  important  at  short  emission  durations  ap- 
plicable to  a  calling  insect,  but  not  at  long  emission  dur- 
ations applicable  to  a  trap.   Second,  near  a  typical 
boundary  surface  the  attractant  concentration  is  much  lower 
than  it  would  be  without  the  boundary,  although  the  concen- 
tration may  remain  at  high  levels  if  the  deposition  velo- 
city  is  extremely  low,  e.g.,  10    cm/sec  in  Fig.  6.   Such  a 
boundary  would  be  considered  a  reflector  rather  than  an 
adsorber.   These  predictions  suggest  that  a  flying  insect 
is  more  likely  to  be  stimulated  by  sex  pheromone  or  other 
attractants  than  an  insect  sitting  or  walking  on  an  adsorp- 
tive  surface.   Observations  by  Visser  (1976)  on  the  host- 
plant  searching  behavior  of  Leptinotarsa  decemlineata  Say 
support  such  a  hypothesis.   A  corollary  hypothesis  is  that 
a  calling  insect  extending  its  pheromone  gland  away  from 
the  surface  on  which  it  is  sitting  has  a  larger  attraction 
space  than  insect  calling  directly  from  the  surface.   The 
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former  calling  behavior  occurs  frequently  in  female  IMM, 
Bombyx  mori  (L. ) ,  and  several  other  insects  (e.g.,  Hammack 
et  al. ,  1976) . 

Consideration  of  plume  shape  effects  in  the  Methods: 
Case  III  section,  led  to  the  hypothesis  that  well-defined 
plumes  from  attractant  sources  of  small  dimensions  are  more 
efficient  in  attracting  insects  than  indistinct  plumes  from 
attractant  sources  of  large  dimensions.   There  are  at  least 
2  experimental  studies  bearing  in  part  on  this  hypothesis. 
Lewis  and  Macaulay  (1976)  found  that  traps  emitting  well- 
defined  smoke  plumes  tended  to  capture  more  insects  than 
those  emitting  indistinct  smoke  plumes.   In  a  related  study 
Macaulay  and  Lewis  (1977)  tested  the  effect  of  source 
dimensions  on  trap  catches  with  inconclusive  results.   How- 
ever ,  there  are  objections  to  the  methodology  of  the  2nd 
study,  in  that  the  large-sized  sources  had  10-fold  to  100- 
fold  greater  emission  rates  than  the  small-sized  sources. 
Because  the  resulting  attraction  spaces  of  the  large-sized 
sources  were  much  larger  than  the  spaces  of  the  small-sized 
sources,  the  attraction  efficiencies  cannot  be  compared  di- 
rectly and  further  study  is  needed  to  confirm  or  deny  the 
hypothesized  plume  shape  effect. 

The  altered-behavior  spaces  predicted  by  the  model  are 
of  theoretical  use,  not  only  for  explaining  decreases  in 
trap  catches  with  increases  in  emission  rate,  but  also  for 
suggesting  the  function  of  the  altered-behavior  threshold. 
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This  is  illustrated  by  the  calculations  in  the  Results: 
Case  II  section,  for  the  altered-behavior  spaces  of  I MM  sex 
pheromone  traps  and  calling  females.   Because  the  entire 
volume  inside  a  100  or  1000  cm  active  space  of  a  trap 
emitting  0.1  ng/sec  was  an  altered-behavior  space,  it  was 
proposed  that  a  trap  emitting  10   ng/sec  might  have  a  grea- 
ter efficiency  in  a  warehouse  than  a  trap  emitting  1  ng/ 
sec.   Such  an  effect  has  not  been  observed  with  IMM,  but 
Vick  et  al.  (1979)  reported  a  similar  effect  in  trapping 
studies  of  Sitotroga  cerealella  (Olivier).   In  a  6.1  x  6.1 
x  2  m  room,  traps  emitting  S.  cerealella  sex  pheromone  at 
0.02  ng/sec  captured  more  males  than  traps  emitting  0.2 
ng/sec. 

The  finding  that  calling  females  of  both  IMM  and  T.  ni 
emit  pheromone  at  rates  sufficient  to  create  small  altered- 
behavior  spaces  suggests  a  function  for  the  altered-behavior 
threshold.   Under  Case  II  dispersal  the  threshold  occurs 
about  6  cm  away  from  an  IMM  female  and  60  cm  from  a  T.  ni 
female.   The  corresponding  distances  can  also  be  calculated 
for  Case  III  dispersal  using  Eq.  II.  6.   In  a  relatively 
calm  airflow  of  50  cm/sec,  the  altered-behavior  threshold 
occurs  about  4  cm  downwind  from  an  IMM  female  and  18  cm 

downwind  from  a  T.  ni  female,  using  K   K   =0.08  cm     , 
—  —  3   y   z  ' 

C  ,  =  1.25  x  10~2  sec/cm3  for  IMM  and  C  ,  =  10~3  sec/cm3 
rd  rd  ' 

for  T.  n_i.   The  existence  of  this  threshold  may  increase 
the  probability  of  a  stimulated  insect  finding  a  calling 
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insect  because  it  alters  the  searching  behavior  from  an 
extensive  search  pattern  to  a  more  intensive  search  pattern 
(Roeiofs,  1978) .   The  limited  range  for  the  altered-behavior 
spaces  of  the  2  insects  supports  this  hypothesis. 

It  remains  to  discuss  hypotheses  explaining  why  search- 
ing insects  tend  to  approach  calling  insects  or  traps  from 
below  in  still  or  nearly  still  air  (Traynier,  1968;  Killiner 
and  Ost,  1971;  Murliss  and  Bettany,  1977) .   The  usual  hypo- 
thesis is  that  pheromone  falls  in  still  air,  but  this  is 
shown  to  be  invalid  in  the  Results:   Case  I  section.   One 
alternative  hypothesis  is  that  the  presence  of  a  reflective 
boundary  below  the  source  could  result  in  a  greater  concen- 
tration of  pheromone  just  above  the  boundary  surface  than  at 
positions  some  distance  above  the  plume  axis.   Inspection  of 
Figs.  3,  4,  and  6  indicates  that  this  could  happen  if  both 
V,  and  t  are  small.   Another  hypothesis  is  that  the  turbu- 
lence level  tends  to  be  proportional  to  height  above  ground; 
thus,  the  plume  may  be  more  well-defined  and  easier  to  fol- 
low near  the  ground  than  at  greater  heights.   However, 
neither  hypothesis  explains  why  an  insect  would  approach  the 
trap  from  below  in  low  airspeeds  but  not  in  high  airspeeds. 
It  is  also  possible  that  the  behavior  is  relatively  indepen- 
dent of  the  actual  pheromone  distribution.   Further  study  is 
necessary  to  resolve  these  questions. 

Considerations  of  predictions  based  on  the  model  leads 
to  the  following  hypotheses  for  optimizing  sex  pheromone 
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trap  design  and  placement:   (1)  Sources  of  small  dimensions 
produce  more  well-defined  plumes  for  longer  distances  than 
sources  of  large  dimensions.   Thus,  trap  catches  could  be 
improved  by  reducing  the  dimensions  of  the  source.   (2)  If 
the  insect  has  an  altered-behavior  threshold,  a  trap  emit- 
ting a  high  rate  and  producing  an  above-threshold  concentra- 
tions could  capture  fewer  insects  than  a  trap  with  a  low 
emission  rate.   Thus,  an  experimental  design  using  a  large 
number  of  traps  with  small  active  spaces  may  be  a  more 
effective  monitor  of  the  pest  population  than  a  design  using 
a  small  number  of  traps  with  large  active  spaces.   A  side 
benefit  of  the  former  design  is  that  it  mitigates  the  ef- 
fects of  temperature  inversion  and  barriers  on  pheromonal 
distribution.   (3)  If  there  are  openings  greater  than  about 
1  cm   in  a  warehouse,  pherbmone  sources  inside  can  attract 
pests  from  outside  the  warehouse.   Warehouse  pest  monitor- 
ing experiments  should  be  designed  to  account  for  this  effect. 


CONCLUSIONS 
A  model  of  insect  attraction  to  odorant  sources  was 
derived,  unique  in  its  treatment  of  adsorptive  boundaries 
obstructing  odorant  flow.   The  generality  and  versatility 
of  the  model  was  demonstrated  by  considering  the  attraction 
of  several  different  insects  to  sex  pheromone  sources  in 
warehouse,  field,  and.  still-air  environments.   Because  the 
important  relationships  of  the  model  are  presented  graphi- 
cally as  well  as  mathematically,  and  because  pheromonal  be- 
havioral  thresholds  generally  range  between  10   and  10 
molecules/cm  ,  the  procedure  for  estimating  attraction  of 
other  insects  to  sex  pheromone  sources  is  manifest.   Calcu- 
lations of  the  model  resulted  in  testable,  experimentally 
supported  hypotheses  regarding  optimal  trap  designs,  and  the 
effects  of  different  parameters  on  the  odorant  concentra- 
tion distribution.   Further  improvement  in  the  precision 
of  the  model  requires  a  better  understanding  of  searching 
behavior,  rather  than  a  more  precise  quantitation  of  odor- 
ant dispersal.   These  considerations  lead  to  the  following 
conclusions:   (1)  The  model  can  be  applied  to  a  wider  range 
of  attraction  problems  than  previous  models,  which  increases 
its  usefulness  for  pest  management  studies  and  theoretical 
analyses  of  the  attraction  process.   (2)  It  offers  new  in- 
sights into  improving  techniques  for  monitoring  insect 
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pest  populations.   (3)  Future  experimental  and  theoretical 
studies  of  insect  attraction  should  be  directed  more  toward 
increasing  the  quantitative' and  qualitative  understanding 
of  searching  behavior  and  less  toward  refining  calculations 
of  the  odorant  concentration  distribution. 


APPENDIX 
DERIVATION  OF  EQUATION  II. 3 

Because  the  mathematics  of  thermal  and  mass  diffusion 
are  equivalent,  Eq.  I I. 3  can  be  derived  from  the  heat- 
transfer  equation  for  a  unit  instantaneous  heat  source 
(Carslaw  and  Jaeger,  1959,  p.  367): 

2  ?    2 

co      (ah'-lT   +   a   g 
1 


27rarr'  a^   +   ah'  (ah'-l)  (A.l) 

n=l  " 


2 

sin(r9    )    sin(r'9    )    exp    (-k6      t) , 
n  n  n 


Where:   T      is  the  temperature,  deg  C; 
a      is  the  boundary  radius,  cm; 
r      is  the  radius  at  which  the  temperature  is 

calculated,  cm; 
r'      is  the  radius  of  the  source,  cm; 
h'      is  the  ratio  H/K,  cm   ; 

H      is  the  coefficient  of  surface  heat  transfer, 

-2     -2     -1     -1 
cal  cm    sec    sec    deg   C; 

K      is  the  thermal  conductivity,  cal  sec-1  cm 

-1 
deg   C; 

2     -1 
<  is  the  thermometric  conductivity,  cm   sec   ; 

9       is  the  nth  positive  root  of  the  equation 

a  cot(ae)  +  (ah'-l)  =  0.  (A. 2) 

GO 
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The  boundary  conditions  specify  that  0  £r<a.  There  is  a 
unit  instantaneous  source  at  r  =  r'.  Heat  transfer  at  a 
is  governed  by  the  equation 

*  =  HT,  (A. 3) 

where:   <j>      is  the  flux,  the  number  of  calories  trans- 
ferred across  a  1  cm"21  surface  in  1  sec 
(Carslaw  and  Jaeger,  1967,  p.  19). 

The  heat  diffusion  equation  can  be  transformed  to  a  mass 
diffusion  equation  by  replacing  T  with  C,  k   with  D,  and  H 
with  V   (Carslaw  and  Jaeger,  1967,  p.  28): 

,    ,     ,,2  2       2 

(ah~l)      +   a      e 
i  n 


2Trarr»  ?      2  (A. 4) 

n=l   a      6      +   ah(ah-l) 
n 

2 

[sin(ro    )sin(r'e    )exp(-De    t) ] , 


where:   h      is  the  ratio,  V  /D,  cm 

d 


To  convert  from  a  diffusion  equation  for  a  unit  instan- 
taneous source  to  the  equation  for  a  source  emitting  at  the 
constant  rate,  Q,  Eq.  A. 4  must  be  integrated  over  the  mea- 
surement period.   It  is  also  convenient  to  use  the  approxi- 
mation: 


6  2 


sin(r'0  ) 

2_   =  1,  (A. 5) 


r'  6n 


which  is  valid  for  r'<<  a.   Accordingly, 

00   (ah-1)2  +  a2  92 


c  =  Af    e    — ; <V  sin  (rV 

n=l  a   6   +  ah (ah-1) 
n 


t        2 
/  exp  [-De    (t-t' ) ]  dt' , 
o        n 


which   gives   Eq.    II.  3    in   the   text: 


(A. 6) 


2  2      2 

»      (ah-1)      +   a      Q 

C   =         Q  E      ^_  (A-7) 

2      2 
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[sin(r6    )/8    ][l-exp(-D      t) ] 
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